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CONJUGATE OUASILINEAR DIRICHLET AND 
NEUMANN PROBLEMS AND A POSTERIORI 
ERROR BOUNDS 


John E. Lavery 

Goddard Space Flight Center 


INTRODUCTION 

Interest in solving quasilinear Dirichlet and Neumann problems of the type 
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arising in such areas as torsion, elasticity, and filtration has been sharpened in the past few 
decades by the development of numerical methods for efficiently solving these problems 
and the possibility of programming these methods on computers. Of great importance for 
the application of these numerical methods are computable a posteriori error bounds for 
the norm of the difference of an approximate solution and the unknown exact solution. 

In this document, the concept of conjugate quasilinear Dirichlet and Neumann problems 
is introduced with the goal of obtaining two-sided bounds for the energies of the solutions 
of the Dirichlet and Neumann problems. These two-sided bounds for the energies are 
subsequently used to obtain a posteriori error bounds for the norm of the difference of 
the approximate and exact solutions.* 

PRELIMINARY REMARKS AND AUXILIARY RESULTS 

All variables and constants will be exclusively real-valued. Lower case c’s, with or without 
numerical subscripts, will denote finite constants. Derivatives will be denoted by literal 
subscripts for functions of two variables, for example, u^^ = 3^u/3y3x (x,y), and by 
numerical subscripts for functions of more than two variables, for example, P 32 (X,Y,x,y) = 
3^p/3Y3x (X,Y,x,y). 

The letter D will denote the rectangular domain (-a, a) x (-b,b) in . The area differential 
in D or in any domain in will be denoted by dA. The boundary D* of D will be para- 
metrized in a counterclockwise direction starting from the lower left corner by a parameter 
s. The four open straight line segments of D* will be denoted by D,*, i = 1 ,2,3,4 (in any order). 

Differential will mean Gateaux differential. Let B and denote Banach spaces. Let B^ 
denote a subspace of B and denote the set {feB: f , f^eB^,}, where is some 

fixed element of B. For operators with differentials, the following generalizations of the 

♦Results presented in this document have been presented in more detail in Chapters I and II of the informal 
publication John E. Lavery, “Solution of a Class of Quasilinear Dirichlet and Neumann Problems by the 
Method of Moments with A Posteriori Error Bounds,” GSFC X-814-73-29, Greenbelt, Maryland, 

January 1973. 
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mean-value theorem and Taylor theorem of classical analysis hold (reference 1 , pp. 37 - 41 ). 
In these lemmas, A denotes an operator from into . 

Lemma 1 . If the operator A possesses a differential A' at every point of Bj , then, for any 
fixed f in Bj , any fixed 77 in B„ , and any fixed e in the conjugate space B| of B^ , there 
exists a constant r, 0 < t < 1 , such that 

(A(f + r?) - A( 0 , e) = (A' (f + t 77, 7?), e), 
where t depends on f, tj, and e. 

Lemma 2 . If the functional A possesses a second differential A" at every point of Bj , and 
A''(f + T77, 77, 77") is continuous as a function of the real variable t for all fixed ? in Bj and 
all fixed 77, t), and 77 in Bq, then, for any fixed f in Bj and any fixed 17 in Bq, there exists a 
constant t, 0 < t < 1 , such that 

A (f + 77) - A (f) = A' (f, 77) + ^ A" (f + T 77, 77, 77). 

A functional A with a second differential will be called convex if there exists a positive 
constant c such that A''(f , 77, 77) > c II 77 Hi . (This definition is more restrictive than usual, 
but is sufficient for our purposes.) 


Lemma 3. Let the Banach space B.be a Hilbert space. Let A be a convex functional on 
B j , and let A satisfy the hypotheses of lemma 2 . Then there exists a unique point f 
in B, at which A has an absolute minimum. The point f is the unique stationary point of 

1 ^ A 

A, that is, the unique point at which A (f ,77) = 0 for all 77 in B^ . 

The Sobolev space W™ (D) will be denoted by H™ (D), with inner product 


(0,0)„m(o) = 


E// 

I i I = 0 D 


^(i) 0(0 (j A. 


L2(D) will denote the Hilbert space of equivalence classes of elements of H^ (D) which 
differ by a constant, with the inner product 


l‘(d) 


•// 


(4>A +^y^y)dA. 


H‘ (D) will denote the Hilbert space which is the closure of the set of infinitely differentiable 
functions with compact support in D in the norm || • 11^1(0) • 

The following lemma is a generalization of a result of classical analysis. (See reference 2 , 
p. 227 , theorem 4 .) 
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Lemma 4. Let 0 and 0 be in (D) and let 


fl 


(0h - 0h ) dA = 0 


for all h in Hj(D). Then there exists a unique in L^fD) such that 

34 > . 34 > 

— =0 and = 0. 

3x 3y 

The standard Friedrichs inequality is stated in lemma 5. 


Lemma 5. For all u in Hj(D), 

llull^, > llull^, . 

[4 min (a^,b^) + 1] ^ 

will denote the Banach space consisting of the traces on D* of the elements of 
(D) (see reference 3). It is well known that C L 2 (D*). For any g e ), 

Hg (D) will denote the set of elements u e (D) such that u|j^. = g (“the trace of u on D’ 
is g”)- Consider the functional M defined for all u in (D) by 



(u^ + up d A. 


This functional satisfies the hypotheses of lemma 3. Therefore, there exists a unique element 
u that minimizes M over (D). The norm II *111/2 Banach space H^/^(D*) will be 

given by 



Let two elements of H^/^(D*) be defined to be equivalent if they differ by a constant. Let 
L^/^(D*) be the Banach space of equivalence classes of defined by this equivalence 

relation, equipped with the norm 


1/2 


(D*) 


Li (D)’ 


where u is the element that minimizes M(u) over Hi (D), g being any member of the equiva- 
lence class g. To each element v in LpD) there corresponds a unique element g = vl^^. in 
Li/^(D*). 


H’1/2(D*) will denote the linear space of linear functionals on H^/^(D’). The value of the 
element g eH“^^^(D*) at geH^/^(D*) will be denoted by (g,g). L^CD*) is a subspace of H'^/^(D*) 


3 


since, to every element g" of LjCD*), there corresponds the functional g of defined 

by 


(g,g) 



gg ds. 


Define (D*) = |g e (g, 1) = o} , Let L*/^(D*) be the linear space of linear 

functionals on L*/^(D’). The value of g e L */^(D*) at g e L*/^(D*) will be denoted by 
(g,g). It can be shown that is isomorphic to L'*/^(D*). For a given g in 

define the functional N« on (D) by 

dA-(g,vlp-)- 

This functional satisfies the hypotheses of lemma 3. Therefore, there exists a unique element 
V that minimizes N; over L:? (D). Define the norm 11*11 on by 

g Z " A / ^ 


Under this norm, is a Banach space. 

The following lemma will permit defining the the derivative of an element of 


Lemma 6. Let g and g" be any two fixed elements of Then the bilinear functional 


(u, h) = 


If. 


(u h -u h)dA 
^ y X X 


has the same value for all (u,h) in (D) X H~(D). 

Let u be the unique element that minimizes M over H* (D). By lemma 6, $(u,h) is a linear 
functional on hlj^. in H^'’^(D*). Thus, the element g' of H'^/^(D’) is well-defined by the 
equality 

(g'.h lo-) = *I>(u,h) 


for all h in H* (D). Moreover, for all h such that hl^. = 1 , 




(g',l)= // [dy(l)^-u^(l)y] dA=0. 


Hence g is in H */^(D”). The element g' of H'*/^(D‘) will be called the derivative of the 
element g of H^''^(D*). Note that g = c implies that u = c and hence that c' = 0. 
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Conversely, the concept of a primitive of an element g of can be defined as 

follows. Let V minimize N| over L^fD). By lemma 3, N|(v,h) = 0 for all h in L‘ (D) or, 
equivalently, for all h in (D). Hence, for all h in Hj(D), 

+ hy)dA = N^(^,h) = 0. 

Lemma 4 then implies that there exists a unique element^u of L^ (D) such that u^ = -Vy and 
Uy Define g = u|p. eL^/^(D*). A primitive of g in H‘^/^(D’)will then be any fixed 
element / g e out of the equivalence class g. Note that any two primitives of g differ 

only by a constant. 

The derivatives and primitives defined here can be multiplied by constants and added in the 
same manner as conventional derivatives and primitives, that is, g' + cg^'= (g + eg")' and 
/g + c/g = Cp + f(g + eg). The following two lemmas also hold. 

Lemma 7. Let g be in H^/^(D*). Then 



f (s') = g + c. 

Lemma 8. Let g be in H *^^(D*). Then 

[(/g) + c]' = g. 

Using lemmas 6 and 7 the following result can be proved. 


Lemma 9. Let g be in H^/^(D*). Then the following two statements are equivalent: iT is in 
(D), and u" is in (D) and is such that 

(iT h - iT h )dA = (g',h L') 

vyx xy^ 



for all h in (D). 

It is of importance that, in all cases of computational interest, derivatives of elements of 
H’/^(D*) and primitives of elements of H'^''^(D*) can be calculated by conventional 
differentiation and integration along D’ rather than by solving the linear Dirichlet and 
Neumann problems of finding the u and v such that M'(u, h) = 0 and Ng (v, h) = 0, that is, 
of minimizing M(u) and Ng(v). These results are stated in the following two lemmas. 

Lemma 10. Let g = g(s) be in|g e C° (D*): gl^- e Wj (Dj), i = 1, 2,3,4 } . By dg/ds, denote 
the element of LjCD*) that equals the distributional derivative (counterclockwise direction) 
of gip. on D.‘, i = 1 ,2,3,4. Then g'= dg/ds. 
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Proof. By lemma 9 and the density of the set of infinitely differentiable functions on D in 
Hi (D). 

Lemmall. Let g = g(s) be in LjfD’) n H‘i/^(D’). Let g = g(s) = g(s) ds. Then 
/g = g + c. 

Proof. By lemmas 7 and 1 0. 

Since the trace on D* of any element of (D) is in L^ (D*), the definition of the Hilbert 
space pHi (D) = {ve H^ (D): ) ds = 0 } with inner product (’jOh^d) made. 

For this space, the following lemma holds. 

Lemma 12. For all v in (D), 

llvll^ > i -llvll\ . 

8max(a2,b2) + l 

To conclude this section, a lemma concerning matrices will be stated. 


Lemma 13. Let 0 (x,y), ^x,y) and 0 (x,y) be elements of Lqo(D), and let there exist positive 
constants c, c, and c" such that 

c < 0 (x, y) < c, 
c < 0 (x, y) < c, 

and 

o' < 0 (x, y) 0 (x, y) - [0' (x, y)] ^ 


for almost all (x,y) in D. Then, for all (X,Y) in , 
(X Y) 

at almost all (x,y) in D. 


! 0 (x, y) 0" (x, y)\ / X 
'0(x,y) 0(x, y)/\ Y 


> — (X^ + Y^) 
2c 


CONJUGATE QUASILINEAR DIRICHLET AND NEUMANN PROBLEMS 

Let there be given functions p(X,Y,x,y) and q(X,Y,x,y) satisfying the following conditions 
(la)-(le); 

p(X, Y, x, y) and q(X, Y, x, y) are defined and finite for all (X, Y, x, y) in X D; \ 
for each fixed (x, y) in D, p(X, Y, x, y) and q(X, Y, x, y) are in C^(R^); /(laj 

for each fixed (X, Y) in R^ , Pj(X, Y, x, y) and q.(X, Y, x, y) are in Loo(D)> i = 1 .2. ) 
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There exist positive constants , m2, M^, and m3 such that for all (X,Y,x,y) in 

r 2 X D 


m, < p(X,Y,x,y)< Mj, 
nij < q (X, Y, X, y) < , 

nij < p (X, Y, X, y) + Pj (X, Y, x. y) X < M^, 
ni2< q(X,Y,x,y) + q2(X,Y,x,y)Y< M^, 


(lb) 


(Ic) 


m3 < [p (X, Y, X, y) + p (X, Y, X, y) X] [q (X, Y, x, y) + q (X, Y, x, y) Y] 

(Id) 

- [P2 (X, Y, X, y) X] [qj (X, Y, X, y) Y] . 

(Condition ( 1 d) guarantees the uniform ellipticity of the Dirichlet and Neumann problems 
and is used to show that the energy functionals for these problems are convex.) For all 
(X,Y,x,y) in X D, 

Pj (X, Y, X, y) X = qj (X, Y, X, y) Y. (jg) 

(This condition is the so-called potentialness condition, which guarantees the existence of 
potentials or energy functionals, the minima of which occur at the solutions of the Dirichlet 
and Neumann problems. See chapter II of reference 1 for general theorems on the existence 
of energy functionals.) Note that conditions (Ic), (Id), and (le) together imply that for all 
(X,Y,x,y) in R^ X D, 


[p (X, Y, X, y) + p (X, Y, X, y) X] [q (X, Y, x, y) + q (X, Y, x, y) Y] 

3 - [p 2 (X,Y,x,y)X] • [qj (X,Y,x,y)Y] 2 (Id) 

Inequalities (Ic) and (Id^) and equality (le) imply that 

Ip^ (X, Y, X, y) XI < 

t I (If) 

Iqj (X,Y,x,y)Yl<yM2-m3 

uniformly for all (X,Y,x,y) in X D. 

Let there be given a function g in H^^^(D*). The Dirichlet problem consists in finding a 
function U e (D) such that 

[p (U^, Uy, x,y) + q (U^, Uy, X, y) h^] dA = 0 (2) 

for all h in (D). (This is the variational formulation of the Dirichlet problem stated in 
the introduction.) 
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Define the energy functional of the Dirichlet problem to be 


£// 


[tp (tu^ , tUy , X, y) + tq (tu^ , tu^ , x, y) ] dA dt 


E(u)=j 

*0 •'“'D 

on Hg (D). Using conditions ( 1 ) it can be shown that the first differential of E is 
E' (u, h) "JJ" tP K> Uy> y) hx + q , x, y) hy] dA 


and that the second differential of E is 


E" (u,h,h) = 


// 


[P K, u X, y) + Pi K, Uy , X, y) uj + [p^ K, Uy, x, y) uj hy 


+ [qj (u^,Uy,x,y) Uy] hy + [q (u^, Uy , x, y) + q^ (u^,Uy,x,y) Uy] hy hy 


> d A, 


where h and h are any elements of Hq(D). The properties of E are used to prove the follow- 
ing two theorems. 

Theorem 1 . There exists a unique solution of the Dirichlet problem. 

Proof. Define the 2-by-2 matrix function P = P(x,y) by 

/p(u^,Uy,x,y) + Pj (u^,Uy,x,y)u^ \ 

P = 

\ qj (Ux>Uy>x,y)Uy q(u^,Uy,x,y) + q2 (u^,Uy,x,y)Uyy. 

From conditions (Ic) and (Id) and lemma 13, one obtains that 


// 


E"(u,h,h)= II (h^,hy)-P- 


d A 


(3) 


5:# 


(h^ + hp d A 


uniformly for all u in (D) and all h in HpD). Lemma 5 further implies that 
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the second equality following from equality (5). By inequality (4), therefore, 

"’a 2 

E(„u)-E(U)> —II u-UII j , (6) 

4Mj [4min(a2,b2)+l] » 

which implies that ^ u ->• U in (D). 

O 

Let there be given a function g in The Neumann problem consists in finding a 

function V e L* (D) such that 

[P Vy, X, y) + q (V^, Vy , x, y) Vyhy] dA - (g, h 1^.) = 0 (7) 
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for all h in L^(D). (This is the variational formulation of the Neumann problem stated in the 
introduction.) 

Define the energy functional of the Neumann problem to be 

F (v) = / 1^1 [tp (tv^ , tVy , X, y) \l + tq (tv^ , tv^ , x, y) ] dAdt - (g, vl^ .) 

D 

on hi (D). Using conditions (1), it can be shown that the first differential of F is 




’ '^x ^ *1 ^'^x ’ '^y > y) ^ ~ (S’ 


and that the second differential of F is 

[p(v^,Vy,x,y)+Pj (v^,Vy,x,y)vJ + [p^ (v^, v^, x,y) vj 


F"(v,h,h) = 


// 


+ hi (v^,Vy,x,y) Vy] hyh^ + [q(v^, v^, x,y) + (v^, v^, x,y) v^] 


d A, 


where h and hare any elements of L^CD). The properties of F are used to prove the follow- 
ing two theorems. 


Theorem 3. There exists a unique solution of the Neumann problem. 


Proof, By the same reasoning as that used to prove inequality (3), one obtains that 


F" (v, h, h) > 


m. 


2M^ 


L^(D) 


( 8 ) 


uniformly for all v and h in L^CD). Hence, F is convex over L^CD). Lemma 3 then implies 
that there exists a unique function V in L^CD) such that 

F' (V, h) = 0 (9) 

for all h in q (D). Let V denote also the unique element in (D) out of the equivalence 
class V which is in (D); Equality (9) then implies that V satisfies equality (7) for all 
h in L^CD) or, equivalently, for all h in (D). Hence, V is the unique solution of the 
Neumann problem. 


Theorem 4. If ^ C LjCD) is a minimizing sequence for F(v), that is, if 

F(„v)-> inf F(v). 
veL^D) 
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as n ->■ oo, then v ^ V in Lj (D). If each element of the sequence is normalized so that 
1 ^ oH^D), then „v - V in (D). 

Proof. By lemma 3, the unique absolute minimum of F(v) occurs at the solution V of the 
Neumann problem, that is, 

inf F(v) = F(V). 
ve Ll (D) 


By lemma 2, for each n > 1, there exists a constant 0<r^ <1, such that 


F („v) - F (V) = F' (V, - V) +- F" (V + 7„ („v - V), „v - V, „v - V) 



(V + T„(„v-V),„v-V, 


-V),' 


the second equality following from equality (9). By inequality (8), therefore. 


F( v)-F(V)> 


m. 


4M, 


ll„v-Vl| 


2 

L^(D)’ 


( 10 ) 


which implies that ->• V in L^(P). If V and the are normalized to be in (D), then 
inequality ( 1 0) and lemma 1 2 imply that 


F(„v)-F(V)> 

4 M 2 [8 max(a^, b^) + 1] 


n 


v-V 


2 

H*(D)’ 


( 11 ) 


which implies that v -> V in ^ (D). 

The concept of conjugate Dirichlet and Neumann problems is based on the following trans- 
formation from to R^. Let functions p(X,Y,x,y) and q(X,Y,x,y) satisfying conditions 
( 1) be given. Let (x,y) be a fixed but arbitrary point in D. Consider the transformation 
(X,Y) -> (Z,W) from R^ into R^ defined by 


Z = q(X, Y. x,y)Y, 
W = -p(X,Y,x,y) X. 


( 12 ) 


This transformation is a special case of the Legendre dual transformation (see reference 4). 
The Jacobian matrix of the transformation is 


qj (X, Y, X, y) Y q (X, Y, x, y) + q^ (X, Y, x, y) Y 


J(X,Y,x,y) = 


-p (X, Y, X, y) - Pi (X, Y, X, y) X - p^ (X, Y, x, y) X 
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a matrix which is nonsingular for all (X,Y) in , since, by condition (Id), 

0< m3 < det J(X,Y,x,y). 

The invertibility of J implies that the range of transformation (12) is open. Using condition 
(lb), it can be shown that the range is also closed. Since the range is nonempty ((0,0) is in 
the range), then, transformation (12) maps R^onto R^ . Let (X,Y) and (X,iO both be 
preimages of some point (Z,W). Then, by lemma 1, 


//z \ /z \ / Y -Y \ \ / ^ /X-X\/Y-Y\ 

y\W/'(w/’(-X + X/ j ^ ~ y \Y-Y/’(-X + X/ 

'p(X,Y,x,y) + Pj (X,Y,x,y)X P2(X,Y,x,y)X 


= (X-X, Y-Y) 


qj (X,Y,x,y) Y 


q(X,Y,x,y) + q,(X,Y,x,y)Y 


X-X 

Y-Y 


>c [(X-}<)2+(Y-Y)2] 


for some positive constant c, waere X = X + t(X - X) and Y = Y + t(Y - Y) for some t in 
(0,1). The inequality above follows from the positive definiteness of the matrix in the 
second-to-last line (see lemma 1 3) and implies that X = X and Y = Y. Hence, transformation 
( 1 2) is not only onto but also one-to-one, that is, it is an isomorphism of R^ onto R^ . This 


permits one to define conjugate functions p(Z,W,x,y) and q(Z,W,x,y) for all (Z,W) in R^ 
by the relations 


p (Z, W, X, y) = 


1 

q(X,Y,x,y) 


q(Z,W, x,y) 


1 

p(X,Y,x, y) 


(13) 


From equalities (12) and (1 3), it is clear that the inverse transformation (Z,W) (X,Y) is 

given by 


X = -(i(Z,W, x,y)W, 
Y = p(Z,W, x,y)Z. 


(14) 


The functions p and q defined by equalities (13) satisfy the same conditions as p and q (the 
constants m. and M. may be different). Indeed, equalities (13) imply that p and q satisfy 
condition (lb) with mj and Mj replaced by m^ = 1/Mj and Mj = 1/mj , respectively. Since 
the mapping (X,Y) ->■ (Z,W) is continuously differentiable with nonsingular Jacobian matrix, 
the inverse mapping is also continuously differentiable. This fact in turn implies that p and 
q satisfy conditions (la), since from equalities (13) one obtains, for example, 
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r 


ax 


3Y 


-qj (X, Y, X, y) — - (X, Y, x, y) — 


Pj (Z, W, X, y) = 

[q(X,Y,x,y)]" 

For fixed (x,y) in D, consider the Jacobian matrix 

- q, (Z, W, X, y)W - q(Z, W, x, y) - q^ (Z, W, x, y)W 


[J(X,Y,x,y)]-» = J(Z,W,x,y) = 


p (Z, W, X, y) + Pj (Z, W, X, y)Z p^ (Z, W, x, y)Z 


of the inverse transformation (14). From condition (Id') for p and q, one obtains condition 
(Id) for p and q: 

nij < [det J (X, Y, x, y)]'* = det J (Z, W, x, y) 


= [p(Z, W, X, y) + Pj (Z, W, X, y)Z] [q(Z, W, x, y) + q^ (Z, W, x, y)W] 
-(p^(Z,W,x,y)Z] [q,(Z,W,x,y)W], 

where nij = l/Mj. Equating the off-diagonal elements of J and J'* and using conditions (Ic) 
and (Id') for p and q yields condition (Ic) for p and q with m 2 and Mj replaced by m 2 = 

= m 2 /(M 2 )^ and M 2 ~ M 2 /m 3 , respectively. Finally, equating the diagonal elements of J and 
J'^ and using condition (le) for p and q imply that p and q satisfy condition ( le). Thus, p 
and q satisfy all of conditions (1) (with different constants). 

The relationship between the energies of the solutions of conjugate Dirichlet and Neumann 
problems is a consequence of the following lemma. 


Lemma 14. For all (X,Y) in , 


p(Z,W,x, y)Z2+q(Z,W, x, y)W^ 



[t p(tZ, tW, X, y) Z^ + 1 q(tZ, tW, x, y)W^] dt 


(15) 



[t p(tX, tY, X, y)X2 + 1 q(tX, tY, X, y)Y^] dt. 


Proof. Fix (x,y) in D. Consider X and Y to be functions of Z and W. Define 0(Z,W) to be 
the left side of equality (15) minus the right side. Using relations (12) and (14), it is easily 
shown that 90/3Z = 90/9W = 0. This, and the fact that 0(0,0) = 0, imply that 0(Z,W) = 0 
for all (Z,W) in , from which equality ( 1 5) follows. 
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Let there be given the Dirichlet problem (2). Define the functions p and q by relations 
(13). Define g in H'^/^(D*) by the equality 


The definition of g' involves solving a linear Dirichlet problem with coefficients p = q = 1 , 
but if g is in |g e C° (D*); gljj. e (D.‘), i = \, 2 , 3 , 4 \ , then, by lemma 10, g' is the con- 
ventional distributional derivaiive of g. The problem of finding V e (D) such that 



x,y)V^h^ + q(V, 


y’ 


x,y)Vyhy] dA-(g,hio.) = o 


(16) 


for all h in (D) is a well-posed Neumann problem with a unique solution. This problem 
will be called the Neumann problem conjugate to the given Dirichlet problem (2). 

Let U be the solution of the Dirichlet problem. Since equality (2) holds for all h in Hq(D), 
lemma 4 implies that there exists a unique v in Lj (D) such that 




Denote also by v the unique element of the equivalence class v which is in (D). Equalities 
(17) imply that 

Uy=p(v^,Vy,x,y)v^. 

Therefore, by lemma 9, 



y) \K Vy , X, y) Vyhy] dA = 




= (g'.h|jj.) 


for all h in (D), that is, v = V, the solution of the Neumann problem (16). Hence, the 
relations 


Vx=q(Ux.Uy.x,y)Uy, 

Vy=-p(U^,Uy,x,y)U^ 


(18) 
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hold between the solution U of the given Dirichlet problem (2) and the solution V of the 
conjugate Neumann problem. 

Let there be given the Neumann problem (16). Define p and q by relations (13). Define 
g in H^^^(D') by the equality 

g = (/g) + c, (19) 

where c is a constant which can be adjusted according to external criteria. The definition of 
f g involves solving a linear Neumann problem with coefficients p = q= 1 , but if g is Lj(D ), 
then, by lemma 1 1 , g coincides with the distributional primitive g(s) = g(s) ds + c of g(s). 
The Dirichlet problem of finding U e Hg(D) such that equality (2) is satisfied, where p, q, 
and g are defined by equalities (13) and (19), will be called the Dirichlet problem conjugate 
to the given Neumann problem ( 1 6). 

Equality ( 1 6) and lemma 4 imply that there exists a unique u in (D) such that 


^=P(V,,V^.x.y)V^. 


(20) 


Denote also by u e (D) the unique element of the equivalence class u which is such that 
/p. ulp. ds = /jj. g ds. Then equalities (20) imply that 



- dA = 


// 


[p (V^ , Vy , X, y) + q (V^ , Vy , x, y) ] dA 


= (g,h|D.) 

for all h in (D). Hence, by equality (19) and lemmas 8 and 9, ulp. = g + c for some con- 
stant c. But since /p. u Ip. ds = /p. g ds (u was chosen so that this condition is satisfied), c 
must be zero, that is, u e (D). Equalities (20) imply that 


^x=9(u,,,Uy.x,y)ay, 

and 

Vy =-p(a,,,{ly,x,y)a^. 


Therefore, by lemma 9, for all h in H‘ (D), 


. J^fP (“x > “y > ^> y) “x^x 9 (u,, , Uy , X, y) Uy hy ] dA = - jTo^yK " ^^x^) 

D *D 

= - ((VId-)',0) 


= 0 , 
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that is, u = U, the solution of the Dirichlet problem (2). Hence, relations ( 1 8) hold be- 
tween the solution V of the given Neumann problem (16) and the solution U of the con- 
jugate Dirichlet problem. 

An important property of the conjugation procedures described above is that the problem 
resulting from two successive conjugations is again the original problem. Specifically, if one 
constructs the Neumann problem conjugate to a given Dirichlet problem and constructs the 
Dirichlet problem conjugate to that Neumann problem, then it is clear from the symmetry 
of relations ( 1 3) and from lemma 7 that one obtains the original Dirichlet problem back again 
(by adjusting the constant c in equality (19) properly). If one constructs the Dirichlet problem 
conjugate to a given Neumann problem and constructs the Neumann problem conjugate to 
that Dirichlet problem, then it is clear from the symmetry of relations ( 1 3) and from lemma 
8 that one obtains the original Neumann problem back again. 

Of greatest importance is the property of conjugate Dirichlet and Neumann problems stated 
in theorem 5. 


Theorem 5. Let the Dirichlet problem (2) and the Neumann problem (16) be conjugate to 
each other and have energy functionals E and F and solutions U and V, respectively. Then 

-F(V) = E(U). (21) 

Remark . Sewell (reference 4) presents a similar result based on the Legendre dual trans- 
formation without discussion of the existence of the potential function of the conjugate 
problem or of the differentiation or integration of the boundary function. 

Proof. Since F'(V,h) = 0 for all h in (D), clearly F'(V,V) = 0 or, equivalently. 




Relations ( 1 8) and lemma 14 (and Fubini’s theorem, which permits interchanging integrations 
over D and over t) then imply that equality (21) holds. 
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A POSTERIORI ERROR BOUNDS 

Let I n=i ^ minimizing sequences for the energy 

functionals of the conjugate Dirichlet and Neumann problems. Then, by theorem 5, the 
inequalities 

- F („ v) < E (U) < E u) (-22a) 

and 

-E(„u)<F(V)<F(^v) (22b) 

hold for all n. Moreover, since Ef^^u) -*■ E(U) and F(jjV) -> F(V), the two-sided error bounds 
in inequalities (22) converge to E(U) and F(V), respectively, as n ^ For the Dirichlet 
solution U, inequalities (6) and (22a) yield the a posteriori error bound 

4M^ [4 min (a^, b^) + 1] 2 

P (.“) + F (n*)] > II n“ - u II „ (23a) 

In certain cases, the a posteriori error bound, 

4M2 2 

[E (^u) + F („v)] > II „u - U II (23b) 

which is derived from inequality (3) in the same manner in which inequality (23a) was derived 
from inequality (4), is more meaningful physically than inequality (23a). For the Neumann 
problem, inequalities (10), (11), and (22b) yield the a posteriori error bounds. 



The author expresses his gratitude to Prof. Ivo Babuska of the University of Maryland for 
numerous fruitful discussions and helpful criticisms. The results presented in this paper were 
obtained while the author was employed by NASA/Goddard Space Flight Center, Greenbelt, 
Maryland. 
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